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Abstract 



In this article, we take the point of view that the .^(1835) be a baryonium 
qj . state and calculate its mass within the framework of the QCD sum rules 

Q \ approach. The numerical value of the mass of the A(1835) is consistent with 



the experimental data. There may be some baryonium component in the 
X(1835) state. 
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§ ■ 1 Introduction 

O 



In 2003, the BES collaboration observed a significant narrow near-threshold en- 
hancement in the proton-antiproton (pp) invariant mass spectrum from the radia- 



tive decay J /if) — > ^pp pp. The enhancement can be fitted with either an S- or 
P-wave Breit-Wigner resonance function. In the case of the S'-wave fitted form, the 
Breit-Wigner mass M = 1859±? ±2 5 MeV and the width T < 30MeU. Recently 
the BES collaboration observed a resonance state X(1835) in the r/7r + 7r _ invari- 
ant mass spectrum in the process J/ip — > 77r + 7r~?/ with the Breit-Wigner mass 
M = (1833.7 ±6.2 ±2. 7) MeV and the width T = (67.7 ±20.3 ±7. 7) MeV [2j. Many 
theoretical works were stimulated to interpret the nature and underlying structures 
of the new particle, there exist many possibilities, for example, the pp bound state, 
the pseudoscalar glueball, non-exotic state, etc [3j HI El E] . 

The X (1835) may be pseudoscalar glueball [3], which can take into account the 
observation of the X(1835) in the r]'ir + 7r~ channel not in the 7r°7r + 7r - channel, while 
the strong coupling to the pp state can be related to the large contribution of the 
gluon axial anomaly to the proton spin. The calculations with lattice QCD and 
QCD sum rules indicate that the pure scalar glueballs lie around (1.5 — 1.7)GeV 
and the pure pseudoscalar glueballs lie around 2.6GeV E], we have to resort to 
special mechanism to pull the mass down to 1.835GeU. 

The X(1835) may also not be exotic state and be the candidate for the second ra- 
dial excitation of the rj' meson, which fulfil the pseudoscalar nonet 7r(1800), i^(1830), 
77(1760), X(1835) [5], such assignment can explain the mass, total decay width, pro- 
duction rate and decay pattern phenomenologically. The decay X(1835) — > i]'7i + ii^ 
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takes place through the emission of a pair of .S-wave it mesons, while the decay 
X(1835) — > r]TT7c has not been observed experimentally yet. Whether or not there 
exists this decay mode is of great importance, further experiments are needed to 
prove or exclude the possibility 

In this article, we take the point of view that the X(1835) be a baryonium with 
the quantum numbers J PC = 0~ + |4j, and calculate its mass in the framework of the 
QCD sum rules approach [5]. The radiative decay of the J /if} is generally believed to 
be glue-rich, which can explain the branching ratio of the decay J /if) — > 77/ (through 
J /if> — > 7GG — > 777') is large while the decay ratio of the J /if) — > 777 (through the 
rj' — rj mixing) is small, about (9.8 ±1.0) x 10~ 4 [10]. The observation of the X(1835) 
in the rf channel not in the 77 channel may be due to the intermediate virtual gluons 
are flavor-neutral and the rj' meson is mainly a SU(3) flavor singlet. The threshold 
2m p = 1876 > 1835 and the width T = 68, the decay X (1835) — > pp takes place 
through the fall apart mechanism with re- arrangement in the color space, while 
suppressed kinematically and the decay occurs only through the tail of the mass 
distribution. 

On the other hand, whether or not there exist some quark configurations which 
can result in the baryonium state is of great importance itself, we explore this 
possibility and propose a special quark configuration, later experimental data can 
confirm or reject this assumption. 

The article is arranged as follows: we derive the QCD sum rules for the mass 
Mx of the X(1835) in section II; in section III, numerical results and discussions; 
section VI is reserved for conclusion. 



2 QCD sum rules for the X(1835) 

In the following, we write down the two-point correlation function II(g 2 ) in the 
framework of the QCD sum rules approach, 

I% 2 ) = z|d 4 xe^(0|T{J 5 (x)J 5 + (0)}|0), (1) 

J 5 (x) = Jki{x)ilbJki{x) , (2) 
Jki{x) = t kmn u^ n {x)C^ a u n {x)^ a di{x) , (3) 
X x = (0|J 5 (0)|X). (4) 

Here the k, I, m, n are color indexes, the C is charge conjunction matrix, the a 
is Lorentz index. The hexaquark states can be classified as "baryonia" if they can 
be described as a single multiquark cluster, or "molecules" if they consist of weakly 
bound NN pairs. We take the pseudoscalar proton-antiproton type interpolating 
current J^{x) to represent the X(1835), if we smear the color indexes, the colored 
constituent Jm{x) has the same structure as the Ioffe current rj(x) which interpolates 
the proton, 

r)(x) = e kmn u^ n (x)C-f a u n (x)-f 5 'j a d k (x), 
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the color indexes k and I in the Jm{x) (in other words, the strong color interactions) 
bind the two constituents as a single multiquark cluster uuduud. If we take the color 
singlet operator r](x) as the basic constituent and choose rj 5 (x) = f/(x)i7 5 ry(x), the 
current 7/5 (x) can interpolate a hexaquark state, whether the compact state or the 
loosely deuteron-like pp bound state, it is difficult to separate the contributions of 
the bound state from the scattering pp state. In this article, we take the X(1835) as a 
baryonium state and choose the current Js(x), although the f]§{x) has non- vanishing 
coupling with the X(1835). 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach [9], we insert a complete series of intermediate states satisfying the 
unitarity principle with the same quantum numbers as the current operator J$(x) 
into the correlation function in Eq.(l) to obtain the hadronic representation. After 
isolating the pole term of the lowest A(1835) state, we obtain the following result, 

In the following, we briefly outline the calculations of the operator product ex- 
pansion in the deep Euclidean space. In order to evaluate the correlation function 
II(g 2 ) at the level of quark-gluon degrees of freedom, we determine the quark prop- 
agator in the presence of the quark and gluon condensates firstljH, 

S ab (x) = (0\T{q a (x)q b (0)}\0) 

i^abX WsGZ i - , * \ $ab(qq) 8 ab (g s qaGq)x 2 

= r TT^WayX + X(J uv ) , 6 

2vr 2 x 4 32n 2 x 2K M " ' 12 192 ' v ; 

where the small masses of the u and d quarks are neglected. Then we substitute 
the quark propagator into the following correlation function to obtain the spectral 
density with the vacuum condensates adding up to dimension- 12, 



n(g 2 ) = 4i (ftxe %q ' x E Hj e kmn e k i vr E Wm , n i r Ti [7 a 7&lpSu>{?) ,> (ar ,> lsrtpS V i(-x)\ 



Tr 



C^'S nH {-x)rCS^J-x)} Tr \c^S ni ,(x)r'CSlJx)] . (7) 



In Eq.(7), we have taken the assumption of the vacuum saturation for the conden- 
sates, the high dimension vacuum condensates are always factorized to lower conden- 
sates with the vacuum saturation in the QCD sum rules, the factorization works well 
in the large iV c limit. It is obvious that such an assumption can not take into account 



2 For the u, d and s quarks, the current masses are small, it is convenient to work in the x- 
representation and adopt the external field method, we follow the routine presented in page-28 and 
page-36 in the last (review) article of Ref.[9] to carry out the operator product expansion. For the 
technical details, one can consult the excellent review "Hadron Properties from QCD Sum Rules", 
L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1. 
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some information in the parameter space, the straight forward calculations with the 
standard operator product expansions can lead to a more general expression for the 
vacuum condensates. We take a simple routine in Eqs.(6-7) to simplify the calcula- 
tion and obtain the result in a special case, it is the common approach to deal with 
the multiquark states with the QCD sum rules, for example, the tetraquark states in 
Ref.[llj. In this article, we take into account the contributions from the quark con- 
densates (qq), gluon condensates ( °* ), mixed condensates {qg s <?Gq), an d neglect 
the contributions from other high dimension condensates which are suppressed by 
large denominators and would not play significant roles. Once the analytical results 
are obtained, then we can take the current-hadron duality below the threshold sq 
and perform the Borel transformation with respect to the variable Q 2 = —q 2 , finally 
we obtain the following sum rule, 

\ 2 x e^l = / dse — , (8) 

Jo 7r 



Imll(s) 9s 7 s 4 (qq) 2 5s(qq) 4 5s 3 (qq)(qg s crGq) 

tt = 2 9 7!7!tt 10 + 2 4 4!4!vr 6 + 3 2 2tt 2 + 2 9 3 2 tt 6 

5s 2 (qg s aGq) 2 s 2 (qq) 2 a s GG s(qq)(qg s aGq) a s GG 

2 10 37r 6 2 6 3% 4 ^ 7T ; 2 8 3% 4 [ 7T ; ' U 

Differentiate the above sum rule with respect to the variable jjt, then eliminate 

Mg 

the quantity Ax, we obtain the QCD sum rule for the mass, 

f° , -T£-ImII(s) p --^Imn(s) 

M 2 X = / dsse — / / dse ^ — . (10) 

Jo 71 Jo ^ 

It is easy to integrate over the variable s, we prefer this formulation for simplicity. 

s 

If we replace the e m b with s n , n = 0, 1, 2, 3, ... , we obtain the finite energy sum 
rule (FESR) [12], 



\* X M£ = l SO dss n1 ^^-, (11) 



M\ = f^^. (12) 

Jo * Jo 11 

The threshold parameter sq is determined by the condition, 

^-M 2 X = 0. (13) 

The FESRs correlate the ground state mass with the continuum threshold sq, and 
separate the ground state from the continuum contributions at the very beginning, 
for some pentaquark currents, there happen exist reasonable stability regions sq [13J. 
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The weight function s n enhances the continuum or the high mass resonances rather 
than the lowest ground state, we must make sure that only the lowest pole terms 
contribute to the FESR below the sq, in some case, a naive stability region sq can 
not guarantee a physically reasonable value of the so [13]. For the hexaquark state, 
the situation is much worse, the stability condition in Eq.(13) can not be satisfied, 
we discard the FESR in this article. 

3 Numerical results and discussions 

The input parameters are taken to be the standard values (uu) = (dd) = (qq) = 
-(0.24±0.01GW) 3 , (qg s aGq) = m 2 (qq),m 2 = (0.8±0.1)GW 2 , (s^) = (0.33GW) 4 
and m u = = 0. In numerical calculation, we observe that the contributions from 
the terms with the gluon condensate (— ) are very small, and neglect the uncer- 
tainty of the gluon condensate. The main contributions to the correlation function in 
Eq.(8) come from the terms with (qq) 2 and (qq)(qg s &Gq), about 85%; the contribu- 
tions from the terms with (qq) 4 and (qg s aGq) 2 are about 35% and 20% respectively, 
and cancel out with each other, the resulting net contributions are less than 15%; 
the contribution comes from the perturbative term is very small, about 1%. We ne- 
glect the contributions from other high dimension condensates which are suppressed 
by large denominators. In the QCD sum rules with the interpolating currents con- 
structed from the multiquark configurations, the main contributions come from the 
terms with the condensates (qq) and (ss) [H], sometimes the mixed condensates 
(qg s <jGq) and (sg s aGs) also play important roles [To]. 

In the following, we discuss the criterion for selecting the threshold parameter 
So and Borel parameter Mb in the QCD sum rules dealing with the multiquark 
states. For the conventional (two-quark) mesons and (three-quark) baryons, the 
hadronic spectral densities are experimentally well known, the separations between 
the ground state and excited states are large enough, the "single-pole + continuum 
states" model works well in representing the phenomenological spectral densities. 
The continuum states can be approximated by the contributions from the asymptotic 
quarks and gluons, and the single-pole dominance condition can be well satisfied, 



here the p pert and p non p stand for the contributions from the perturbative and non- 
perturbative part of the spectral density respectively. From the condition in Eq.(14), 
we can obtain the maximal value of the Borel parameter M^ ax , exceed this value, 
the single-pole dominance will be spoiled. On the other hand, the Borel parameter 
must be chosen large enough to warrant the convergence of the operator product 
expansion and the contributions from the high dimension vacuum condensates which 
are known poorly are of minor importance, the minimal value of the Borel parameter 
M™ n can be determined. 




(14) 
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For the conventional mesons and baryons, the Borel window Mg ax — M™ n is 
rather large and the reliable QCD sum rules can be obtained. However, for the 
multiquark states i.e. tetraquark states, pentaquark states, hexaquark states, etc, 
the spectral densities p ~ s n with n is larger than the ones for the conventional 

hadrons, the integral J °° s n e m b ds converges more slowly [13]. If one do not want to 
release the condition in Eq.(14), we have to either postpone the threshold parameter 
So to very large values or choose very small values of the Borel parameter M^ ax . 
With large values of the threshold parameter s , for example, s 3> M 2 r , here the 
gr stands for the ground state, the contributions from the high resonance states and 
continuum states are included in, we can not use the single-pole (or ground state) 
approximation for the spectral densities; on the other hand, with very small values of 
the Borel parameter Mg ax , the operator product expansion is broken down, and the 
Borel window M^ ax — M™ n shrinks to zero or negative values. We should resort to 
the "multi-pole + continuum states" to approximate the phenomenological spectral 
densities. The onset of the continuum states is not abrupt, the ground state, the first 
excited state, the second excited state, etc, the continuum states appear sequentially; 
the excited states may be loose bound states and have large widths. The threshold 
parameter sq is postponed to large value, at that energy scale, the spectral densities 
can be well approximated by the contributions from the asymptotic quarks and 
gluons, and of minor importance for the sum rules. 

The present experimental knowledge about the phenomenological hadronic spec- 
tral densities of the multiquark states is rather vague, even the existence of the 
multiquark states is not confirmed with confidence, and no knowledge about either 
there are high resonances or not. 

In this article, the following criteria are taken. We choose the suitable val- 
ues of the Borel parameter Mb, on the one hand the minimal values Mj^ m are 
large enough to warrant the convergence of the operator product expansion, for 
M™ n > \/3.5GeV, the dominating contributions come from the terms with (qq) 2 
and {qq)(qg s cGq), about 85%; on the other hand the maximal values M^ ax are 
small enough to suppress the contributions from the high excited states and con- 
tinuum states, we choose the naive analysis e~ s °^ M B ax ) < eT 1 . For the hadronic 
spectral density, the more phenomenological analysis is preferred, we approximate 
the spectral density with the contribution from the single-pole term, the threshold 
parameter s is taken slightly above the ground state mass ( s/so > M gr + -^-) to 
subtract the contributions from the excited states and continuum states. In this 
article, the threshold parameter sq is taken to be ^/sq = (2.1 — 2.3)GeV > 1.9GeV. 
It is reasonable for the Breit-Wigner mass M = (1833.7 ± 6.2 ± 2. 7)MeV and width 
T = (67.7 ± 20.3 ± 7.7)MeV. The values of the from the sum rules in Eq.(8) 
increase quickly with ^ l[ /s^ > 2.3GeV, which are shown in Fig.l, it may serve as in- 
dication of the onset of the high resonances and continuum states. From the Fig.l, 
we can see that the Borel parameter can be chosen to be M| = (3.5 - 5.5)GeV 2 , 
M% ax < ^ < 2.3GeV. 
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Figure 1: The M x and with the Borel parameter M\ for (gg) = — (0.24GeV) 3 
and = 0.8GeV 2 . 



Finally, we obtain the value of the mass of the X(1835), 

M x = (1.9±0.1)GeV . (15) 

The numerical result is compatible with experimental data, one may reject taking 
the value from the more phenomenological analysis as quantitatively reliable, the 
result is qualitative at least. The systematic studies with the random instanton 
liquid model indicate that the masses of the diquarks are = (420 ± 30)MeV, 
m v = m A = (940 ± 20)MeV, m T = (570 ± 20)MeV^, here the S, V } A and 
T stand for the scalar, vector, axial-vector and tensor diquarks respectively. In 
this article, the chosen quark configuration has two diquark constituents, a diquark 
(^abcU^Cj^Uc) and an antidiquark (eabcUbC^^u^), it is not surprise that the energy 
scale set by the diquarks is about l.QGeV. Different quark configurations can result 
in different masses for the hadrons, for example, the meson-meson type interpolating 
currents indicate the masses of the tetraquark states are less than lGeV |17j . 



4 Conclusion 

In this article, we take the point of view that the X(1835) be a baryonium state 
and calculate its mass within the framework of the QCD sum rule approach. The 
numerical value of the mass of the X(1835) is consistent with the experimental data. 
There may be some baryonium component in the X(1835) state. 
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